Quantitative Analysis of Single Nucleotide Polymorphisms within Copy Number Variation

Background Single nucleotide polymorphisms (SNPs) have been used extensively in genetics and epidemiology studies. Traditionally, SNPs that did not pass the Hardy-Weinberg equilibrium (HWE) test were excluded from these analyses. Many investigators have addressed possible causes for departure from HWE, including genotyping errors, population admixture and segmental duplication. Recent large-scale surveys have revealed abundant structural variations in the human genome, including copy number variations (CNVs). This suggests that a significant number of SNPs must be within these regions, which may cause deviation from HWE. Results We performed a Bayesian analysis on the potential effect of copy number variation, segmental duplication and genotyping errors on the behavior of SNPs. Our results suggest that copy number variation is a major factor of HWE violation for SNPs with a small minor allele frequency, when the sample size is large and the genotyping error rate is 0~1%. Conclusions Our study provides the posterior probability that a SNP falls in a CNV or a segmental duplication, given the observed allele frequency of the SNP, sample size and the significance level of HWE testing.

Genome-wide scan of copy number variation in late-onset Alzheimer's disease.

Alzheimer's disease is a complex and progressive neurodegenerative disease leading to loss of memory, cognitive impairment, and ultimately death. To date, six large-scale genome-wide association studies have been conducted to identify SNPs that influence disease predisposition. These studies have confirmed the well-known APOE epsilon4 risk allele, identified a novel variant that influences disease risk within the APOE epsilon4 population, found a SNP that modifies the age of disease onset, as well as reported the first sex-linked susceptibility variant. Here we report a genome-wide scan of Alzheimer's disease in a set of 331 cases and 368 controls, extending analyses for the first time to include assessments of copy number variation. In this analysis, no new SNPs show genome-wide significance. We also screened for effects of copy number variation, and while nothing was significant, a duplication in CHRNA7 appears interesting enough to warrant further investigation.

Copy number variation and Huntington’s disease

Tese de doutoramento, Ciências Biomédicas (Neurociências), Universidade de Lisboa, Faculdade de Medicina, 2014

Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

Copy number variation and transposable elements feature in recent, ongoing adaptation at the Cyp6g1 locus

The increased transcription of the Cyp6g1 gene of Drosophila melanogaster, and consequent resistance to insecticides such as DDT, is a widely cited example of adaptation mediated by cis-regulatory change. A fragment of an Accord transposable element inserted upstream of the Cyp6g1 gene is causally associated with resistance and has spread to high frequencies in populations around the world since the 1940s. Here we report the existence of a natural allelic series at this locus of D. melanogaster, involving copy number variation of Cyp6g1, and two additional transposable element insertions (a P and an HMS-Beagle). We provide evidence that this genetic variation underpins phenotypic variation, as the more derived the allele, the greater the level of DDT resistance. Tracking the spatial and temporal patterns of allele frequency changes indicates that the multiple steps of the allelic series are adaptive. Further, a DDT association study shows that the most resistant allele, Cyp6g1-[BP], is greatly enriched in the top 5% of the phenotypic distribution and accounts for ~16% of the underlying phenotypic variation in resistance to DDT. In contrast, copy number variation for another candidate resistance gene, Cyp12d1, is not associated with resistance. Thus the Cyp6g1 locus is a major contributor to DDT resistance in field populations, and evolution at this locus features multiple adaptive steps occurring in rapid succession.

Copy number variation of individual cattle genomes using next-generation sequencing

Copy number variations (CNVs) affect a wide range of phenotypic traits; however, CNVs in or near segmental duplication regions are often intractable. Using a read depth approach based on next-generation sequencing, we examined genome-wide copy number differences among five taurine (three Angus, one Holstein, and one Hereford) and one indicine (Nelore) cattle. Within mapped chromosomal sequence, we identified 1265 CNV regions comprising similar to 55.6-Mbp sequence-476 of which (similar to 38%) have not previously been reported. We validated this sequence-based CNV call set with array comparative genomic hybridization (aCGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH), achieving a validation rate of 82% and a false positive rate of 8%. We further estimated absolute copy numbers for genomic segments and annotated genes in each individual. Surveys of the top 25 most variable genes revealed that the Nelore individual had the lowest copy numbers in 13 cases (similar to 52%, chi(2) test; P-value <0.05). In contrast, genes related to pathogen- and parasite-resistance, such as CATHL4 and ULBP17, were highly duplicated in the Nelore individual relative to the taurine cattle, while genes involved in lipid transport and metabolism, including APOL3 and FABP2, were highly duplicated in the beef breeds. These CNV regions also harbor genes like BPIFA2A (BSP30A) and WC1, suggesting that some CNVs may be associated with breed-specific differences in adaptation, health, and production traits. By providing the first individualized cattle CNV and segmental duplication maps and genome-wide gene copy number estimates, we enable future CNV studies into highly duplicated regions in the cattle genome.

Germline DNA copy number variation in familial and early-onset breast cancer

Introduction: Genetic factors predisposing individuals to cancer remain elusive in the majority of patients with a familial or clinical history suggestive of hereditary breast cancer. Germline DNA copy number variation (CNV) has recently been implicated in predisposition to cancers such as neuroblastomas as well as prostate and colorectal cancer. We evaluated the role of germline CNVs in breast cancer susceptibility, in particular those with low population frequencies (rare CNVs), which are more likely to cause disease." Methods: Using whole-genome comparative genomic hybridization on microarrays, we screened a cohort of women fulfilling criteria for hereditary breast cancer who did not carry BRCA1/BRCA2 mutations. Results: The median numbers of total and rare CNVs per genome were not different between controls and patients. A total of 26 rare germline CNVs were identified in 68 cancer patients, however, a proportion that was significantly different (P = 0.0311) from the control group (23 rare CNVs in 100 individuals). Several of the genes affected by CNV in patients and controls had already been implicated in cancer. Conclusions: This study is the first to explore the contribution of germline CNVs to BRCA1/2-negative familial and early-onset breast cancer. The data suggest that rare CNVs may contribute to cancer predisposition in this small cohort of patients, and this trend needs to be confirmed in larger population samples.

MSX2 copy number increase and craniosynostosis: copy number variation detected by array comparative genomic hybridization

Fleury Group

Single-Nucleotide Polymorphism and Copy Number Variation of the Multidrug Resistance-1 Locus of Plasmodium vivax: Local and Global Patterns

Emerging resistance to chloroquine (CQ) poses a major challenge for Plasmodium vivax malaria control, and nucleotide substitutions and copy number variation in the P. vivax multidrug resistance 1 (pvmdr-1) locus, which encodes a digestive vacuole membrane transporter, may modulate this phenotype. We describe patterns of genetic variation in pvmdr-1 alleles from Acre and Amazonas in northwestern Brazil, and compare then with those reported in other malaria-endemic regions. The pvmdr-1 mutation Y976F, which is associated with CQ resistance in Southeast Asia and Oceania, remains rare in northwestern Brazil (1.8%) and its prevalence mirrors that of CO resistance worldwide. Gene amplification of pvmdr-1, which is associated with mefloquine resistance but increased susceptibility to CO, remains relatively rare in northwestern Brazil (0.9%) and globally (< 4%), but became common (> 10%) in Tak Province, Thailand, possibly because of drug-mediated selection. The global database we have assembled provides a baseline for further studies of genetic variation in pvmdr-1 and drug resistance in P. vivax malaria.

Copy number variation in Williams-Beuren syndrome: suitable diagnostic strategy for developing countries

Abstract Background Williams-Beuren syndrome (WBS; OMIM 194050) is caused by a hemizygous contiguous gene microdeletion at 7q11.23. Supravalvular aortic stenosis (SVAS), mental retardation, and overfriendliness comprise typical symptoms of WBS. Although fluorescence in situ hybridization (FISH) is considered the gold standard technique, the microsatellite DNA markers and multiplex ligation-dependent probe amplification (MLPA) could be used for to confirm the diagnosis of WBS. Results We have evaluated a total cohort of 88 patients with a suspicion clinical diagnosis of WBS using a collection of five markers (D7S1870, D7S489, D7S613, D7S2476, and D7S489_A) and a commercial MLPA kit (P029). The microdeletion was present in 64 (72.7%) patients and absent in 24 (27.3%) patients. The parental origin of deletion was maternal in 36 of 64 patients (56.3%) paternal in 28 of 64 patients (43.7%). The deletion size was 1.55 Mb in 57 of 64 patients (89.1%) and 1.84 Mb in 7 of 64 patients (10.9%). The results were concordant using both techniques, except for four patients whose microsatellite markers were uninformative. There were no clinical differences in relation to either the size or parental origin of the deletion. Conclusion MLPA was considered a faster and more economical method in a single assay, whereas the microsatellite markers could determine both the size and parental origin of the deletion in WBS. The microsatellite marker and MLPA techniques are effective in deletion detection in WBS, and both methods provide a useful diagnostic strategy mainly for developing countries.

MYC, FBXW7 and TP53 copy number variation and expression in Gastric Cancer

Abstract Background MYC deregulation is a common event in gastric carcinogenesis, usually as a consequence of gene amplification, chromosomal translocations, or posttranslational mechanisms. FBXW7 is a p53-controlled tumor-suppressor that plays a role in the regulation of cell cycle exit and reentry via MYC degradation. Methods We evaluated MYC, FBXW7, and TP53 copy number, mRNA levels, and protein expression in gastric cancer and paired non-neoplastic specimens from 33 patients and also in gastric adenocarcinoma cell lines. We also determined the invasion potential of the gastric cancer cell lines. Results MYC amplification was observed in 51.5% of gastric tumor samples. Deletion of one copy of FBXW7 and TP53 was observed in 45.5% and 21.2% of gastric tumors, respectively. MYC mRNA expression was significantly higher in tumors than in non-neoplastic samples. FBXW7 and TP53 mRNA expression was markedly lower in tumors than in paired non-neoplastic specimens. Moreover, deregulated MYC and FBXW7 mRNA expression was associated with the presence of lymph node metastasis and tumor stage III-IV. Additionally, MYC immunostaining was more frequently observed in intestinal-type than diffuse-type gastric cancers and was associated with MYC mRNA expression. In vitro studies showed that increased MYC and reduced FBXW7 expression is associated with a more invasive phenotype in gastric cancer cell lines. This result encouraged us to investigate the activity of the gelatinases MMP-2 and MMP-9 in both cell lines. Both gelatinases are synthesized predominantly by stromal cells rather than cancer cells, and it has been proposed that both contribute to cancer progression. We observed a significant increase in MMP-9 activity in ACP02 compared with ACP03 cells. These results confirmed that ACP02 cells have greater invasion capability than ACP03 cells. Conclusion In conclusion, FBXW7 and MYC mRNA may play a role in aggressive biologic behavior of gastric cancer cells and may be a useful indicator of poor prognosis. Furthermore, MYC is a candidate target for new therapies against gastric cancer.

MYC, FBXW7 and TP53 copy number variation and expression in gastric cancer

BACKGROUND: MYC deregulation is a common event in gastric carcinogenesis, usually as a consequence of gene amplification, chromosomal translocations, or posttranslational mechanisms. FBXW7 is a p53-controlled tumor-suppressor that plays a role in the regulation of cell cycle exit and reentry via MYC degradation. METHODS: We evaluated MYC, FBXW7, and TP53 copy number, mRNA levels, and protein expression in gastric cancer and paired non-neoplastic specimens from 33 patients and also in gastric adenocarcinoma cell lines. We also determined the invasion potential of the gastric cancer cell lines. RESULTS: MYC amplification was observed in 51.5% of gastric tumor samples. Deletion of one copy of FBXW7 and TP53 was observed in 45.5% and 21.2% of gastric tumors, respectively. MYC mRNA expression was significantly higher in tumors than in non-neoplastic samples. FBXW7 and TP53 mRNA expression was markedly lower in tumors than in paired non-neoplastic specimens. Moreover, deregulated MYC and FBXW7 mRNA expression was associated with the presence of lymph node metastasis and tumor stage III-IV. Additionally, MYC immunostaining was more frequently observed in intestinal-type than diffuse-type gastric cancers and was associated with MYC mRNA expression. In vitro studies showed that increased MYC and reduced FBXW7 expression is associated with a more invasive phenotype in gastric cancer cell lines. This result encouraged us to investigate the activity of the gelatinases MMP-2 and MMP-9 in both cell lines. Both gelatinases are synthesized predominantly by stromal cells rather than cancer cells, and it has been proposed that both contribute to cancer progression. We observed a significant increase in MMP-9 activity in ACP02 compared with ACP03 cells. These results confirmed that ACP02 cells have greater invasion capability than ACP03 cells. CONCLUSION: In conclusion, FBXW7 and MYC mRNA may play a role in aggressive biologic behavior of gastric cancer cells and may be a useful indicator of poor prognosis. Furthermore, MYC is a candidate target for new therapies against gastric cancer.

Copy number variation of the Beta-defensin genes in europeans: no supporting evidence for association with lung function, chronic obstructive pulmonary disease or asthma

Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02-1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72-1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.

Extensive Copy-Number Variation of Young Genes across Stickleback Populations

Duplicate genes emerge as copy-number variations (CNVs) at the population level, and remain copy-number polymorphic until they are fixed or lost. The successful establishment of such structural polymorphisms in the genome plays an important role in evolution by promoting genetic diversity, complexity and innovation. To characterize the early evolutionary stages of duplicate genes and their potential adaptive benefits, we combine comparative genomics with population genomics analyses to evaluate the distribution and impact of CNVs across natural populations of an eco-genomic model, the three-spined stickleback. With whole genome sequences of 66 individuals from populations inhabiting three distinct habitats, we find that CNVs generally occur at low frequencies and are often only found in one of the 11 populations surveyed. A subset of CNVs, however, displays copy-number differentiation between populations, showing elevated within-population frequencies consistent with local adaptation. By comparing teleost genomes to identify lineage-specific genes and duplications in sticklebacks, we highlight rampant gene content differences among individuals in which over 30% of young duplicate genes are CNVs. These CNV genes are evolving rapidly at the molecular level and are enriched with functional categories associated with environmental interactions, depicting the dynamic early copy-number polymorphic stage of genes during population differentiation.